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Abstract 



We determine the two-loop 'time-like' Altarelli-Parisi splitting functions, appearing in 
the next-to-leading order Q 2 -evolution equations for fragmentation functions, via ana- 
lytic continuation of the corresponding 'space-like' splitting functions for the evolution of 
parton distributions. We do this for the case of unpolarized fragmentation functions and 
- for the first time - also for the functions describing the fragmentation of a longitudi- 
nally polarized parton into a longitudinally polarized spin- 1/2 hadron such as a A baryon. 
Our calculation is based on the method proposed and employed by Curci, Furmanski and 
Petronzio in the unpolarized case in which we confirm their results. 



1 Introduction 



The spin structure of longitudinally polarized nucleons has been investigated in a num- 
ber of experiments |l) in recent years by scattering polarized highly virtual space-like 
(q 2 = — Q 2 < 0) photons off polarized targets, which provides access to the spin-dependent 
'space-like' parton distributions of the nucleon. In contrast to this, nothing is known as yet 
about the corresponding polarized 'time-like' parton densities, i.e., the functions describ- 
ing the fragmentation of a longitudinally polarized quark or gluon into a longitudinally 
polarized (spin-1/2) hadron. In analogy with the space-like case these are defined by 

AD h f (z,Q 2 ) = D h f §(z,Q 2 ) - D^-^Q 2 ) , (1) 

where D^^(z, Q 2 ) (D^ + ^(z,Q 2 )) is the probability for finding a hadron h with positive 
(negative) helicity in a parton / with positive helicity at a mass scale Q, carrying a 
fraction z of the parent parton's momentum. We note that taking the sum instead of the 
difference in Eq. ([!]) one obtains the corresponding unpolarized fragmentation function. 

Spin-dependent fragmentation functions appear equally interesting as the space-like 
distributions since they obviously contain information on how the spin of a fragmenting 
parton is transmitted to that of the produced hadron. The most likely candidate for 
a measurement of polarized fragmentation functions is the A baryon since its dominant 
decay A — > p7t~ is parity-violating and enables the determination of the A's polarization 
0. In H a strategy was proposed for extracting the A/}j (/ = q,q) in single-particle 
inclusive e + e~ annihilation (SIA) e + e~ — > AX. If the energy is far below the Z resonance, 
one longitudinally polarized beam is required in order to fix the polarization of the out- 
going (anti) quark that fragments into the A and to obtain a non- vanishing twist-two spin 
asymmetry. At higher energies, no beam polarization is needed since the parity-violating 
coupling qqZ automatically generates a net polarization of the quarks. Apart from SIA, 
which plays the same fundamental role for the determination of fragmentation functions 
as deep-inelastic scattering (DIS) does for that of space-like parton distributions, the pos- 
sibility of extracting the AZ}j in semi-inclusive DIS (SIDIS) in the current fragmentation 
region, ep — > eAX, has also been studied theoretically recently |4], [5|. Here, either a 
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longitudinally polarized lepton beam or a polarized nucleon target would be required. On 
the experimental side, ALEPH has reported a first measurement of the A polarization in 
Z decays recently ||, and a measurement of the AD A in SIDIS appears possible for the 
HERMES experiment and is planned by the COMPASS collaboration g. 

The polarized cross section for, sayQ, the process e + e~ — > AX (the arrows denoting 
longitudinal polarization) can be written in the factorized form as a sum over convolutions 
of polarized hard subprocess cross sections with the process-independent fragmentation 
functions ADf of (|), 

dAa e+e ~^ AX (s,y,z) _ da(X e = +, A A = +, s,y, z) da(X e = +, A A = -, s,y, z 



dydz dydz dydz 



(2) 



where A e , Aa denote the helicities of the polarized electron and the A, and s = 2p e + -p e -, 
z = 2pA • q/Q 2 with the momentum q of the time-like (q 2 = Q 2 > 0) intermediate 7 or Z. 
The variable y is defined by y = p^-p e - /pa'Q an d is related to the cms scattering angle 9 of 
the produced A with respect to the incoming electron via y = (1 + cos 0)/2. The polarized 
subprocess cross sections dAa e+e ^^ x (s,y,C,)/dyd^ are defined in complete analogy with 
(0), and £ is the partonic counterpart of z, £ = 2pj ■ q/Q 2 . Again the corresponding 
expression for the unpolarized cross section is obtained by taking the sum instead of the 
difference in Eq. (0) and omitting all A's. Unlike the fragmentation functions, the hard 
subprocess cross sections are calculable in perturbative QCD. QCD can, however, predict 
the Q 2 dependence of the fragmentation functions via the Altarelli-Parisi equations M, 
once a suitable non-perturbative hadronic input for the evolution has been found. In the 
leading order (LO), there is only one subprocess, namely e + e~ — > qq via 7 or Z exchange, 
and the fragmentation functions evolve according to the LO polarized (time-like) Altarelli- 
Parisi equations. 

It is the main purpose of this paper to set up the complete next-to-leading order 
(NLO) framework of QCD for single-inclusive annihilation into a polarized hadron. This 

1 Since the A appears the most realistic candidate for measurements of spin-dependent fragmentation 
functions we write down all formulae below for this specific case. They apply, of course, equally well to 
any other final state (spin-1/2) hadron whose longitudinal polarization can be determined experimentally. 
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task first of all involves calculating the 0(a s ) corrections to the LO hard subprocess cross 
section, including the calculation of subprocesses that first appear at NLO. This was 



recently achieved in |L0|, pj| . However, knowledge of the underlying hard subprocesses to 
NLO accuracy cannot be the full story as becomes immediately obvious from the well- 
known fact that the corresponding corrections are factorization scheme dependent, i.e., 
depend on the convention adopted when subtracting the collinear singularities appearing 
in the calculation. For a fully consistent NLO calculation one also needs to perform the 
evolution of the polarized fragmentation functions in NLO, which requires knowledge of 
the polarized NLO evolution kernels in the (time-like) Altarelli-Parisi equations. Only 
when both types of NLO corrections, those to the subprocess cross sections and to the 
evolution kernels, are known does the NLO framework become complete and consistent, 
the factorization scheme dependencies cancelling out to the order considered whenever a 
physical cross section is calculated. This situation is of course completely the same as in 
the more familiar space-like case of, e.g., DIS structure functions. 

As will be discussed below, it is possible to derive the polarized NLO time-like evo- 
lution kernels by analytic continuation of their space-like counterparts which have been 
calculated recently O, 13 1. The procedure for doing this has first been worked out for 



the unpolarized non-singlet case in JT3J and has also been used for the unpolarized singlet 



sector in flSfl . We will pursue this method. The results we obtain refer to the MS scheme 
and need to be combined with NLO corrections to the hard subprocess cross sections in 
the same scheme, as recently presented for SIA and SIDIS in |TT|. Since the procedure 
of analytic continuation can also be applied to the hard subprocess cross sections we will 
provide a check on the results of for SIA. 



In view of the present lack of any experimental information on the AD^ one could 
argue that it is somewhat premature to set up the full NLO framework for their evolution 
and the processes in which they appear. On the other hand, it seems likely that data will 
become available in the future. Furthermore, the transition from the space-like to the 
time-like region in the polarized case appears interesting in itself: In LO the space-like 
and time-like Altarelli-Parisi evolution kernels are related to each other via an analytic 



3 



continuation rule (ACR) |I7J and also via the so-called Gribov-Lipatov relation (GLR) 
T8| . In the unpolarized case these relations were shown to be broken beyond leading order 



in the MS scheme |L5], [16], and a similar feature is thus expected for the spin-dependent 
case. The NLO effects also appear interesting from a more physical point of view. For 
instance, one would expect [^IJ that to a first approximation polarized-A production in 
SIA proceeds just via strange quark fragmentation s — > A, i.e., is essentially sensitive 
to AD^. NLO evolution on the other hand automatically generates non-vanishing non- 
strange fragmentation functions due to the existence of flavor non-diagonal quark-to-quark 
splitting functions. Also, the possibly important JTDJ role played by gluons is appreciated 
when going beyond the leading order. 

The remainder of this paper is organized as follows. In section 2 we set the general 
framework for our calculations and briefly discuss the LO results. In sections 3 and 4 
we present in some detail the determination of the NLO corrections for the unpolarized 
time-like situation by analytic continuation of their space-like counterparts. Even though 
neither the method of analytic continuation nor the final result of the calculation are new, 
the full calculation itself has never been documented before, and we also provide new 
insight in the breakdown of the ACR beyond LO. Furthermore, our findings in sections 
3,4 are crucial for dealing with the polarized case, which is then done in the subsequent 
section. In section 6 we study an interesting supersymmetric relation obeyed by the NLO 
unpolarized and polarized time-like splitting functions. Section 7 briefly summarizes our 
work. 



2 General framework and LO results 



Let us first set the notation by collecting all ingredients for a NLO treatment of the cross 
section in Eq. We begin by dealing with the hard subprocess cross sections. In 

analogy with the familiar space-like g\ = g[ S \x, Q 2 ) (where x = Q 2 /2p ■ q < 1) we define 
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a time-like structure function g± (z, Q 2 ) and write Eq. (^|) as0 



dAa" 



*AX, 



Qua 



(2y-l)g?\z,Q 2 ). 



(4) 



dzdy Q 2 

To facilitate the further discussion, we adopt a combined treatment of the space-like and 
time-like situations and introduce the structure function G[ U \£,,Q 2 ), where the index 
U stands for either 'space-like' (U = S, Q[ S ' = 2g[ s \ £ = x), or 'time-like' (U = T, 



(T) 



(T) c 

9i , € 



and the parton distributions Af^ u \^,Q 2 ) (/ = q,q,g), where 



A/( s ) = A/ (with the usual polarized hadronic parton densities A/) and Af^ = AD^. 
In terms of the A/^ we can write (?j i7 ' ) to NLO as 

Sr(£,Q 2 ) =Ee, 2 I [A^) + A^] ® ACf)+^A^)® ACf)|(e,g 2 ) , (5) 

where the sum runs over the n/ active quark flavors, rjs = ^/ n f, Vt — 2 and £g> denotes 
the usual convolution. The hard subprocess cross sections AC^ , AC^O are taken to have 
the perturbative expansion 

,(U) > ~ ^ * ~ (rn cm / a 



Acr(e,< 

where ACf >-(°>(£) = 5(1 - £), ACfO-(°)(£) = 0. 



Acf ),(o) (0 + ^ Acf ),(D (0 



(6) 



To determine the Q 2 evolution of the space-like and time-like parton densities Af^ 
in Eq. (^) it is as usual convenient to decompose them into flavor singlet and non-singlet 
pieces by introducing the densities Ag±^ and the vector 

^-(^). 

where 

Agf^Ag^iAg^, AE (C7) = ^(Ag (c/) + Ag (c/) ) . (8) 

One then has the following non-singlet evolution equations (q, q being two different fla- 
vors) : 



d 



(Agf-AgH^Q 



dlnQ s 



,„, + ® (Agr - Agf >) 



d 



dlnQ 1 



(9) 
(10) 



2 For simplicity we restrict our considerations to pure photon exchange in the process e + e — » gg. 
Exchange of Z° and "fZ° interference modify the angular dependence of the longitudinally polarized 



cross section and thus add new structure functions to its expression |10 



The two evolution kernels APgF±(£, a s (Q 2 )) start to become different beyond LO as a 
result of the presence of transitions between quarks and antiquarks. The singlet evolution 
equation reads 

d 



AvW(Z,Q 2 )= ApW®AtfW (£,Q 2 ). 



d\nQ 2 

We write the singlet evolution matrices for the space-like and time-like cases as 

APi T ) 2n f AP<P 
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AP( S ) AP( S ) 

99 99 

AP(S) ap(s) 

99 99 



The gg-entries in (1121) are expressed as 



ApCO 



99 



99 



h^ap(j) APir) 



(12) 



2n 



99 



99 



^ r qq ~ ^ r qq,+ + ^^qq^PS ■ 



(13) 



APqq p S which vanishes in LO is called the 'pure singlet' splitting function since it only ap- 



pears in the singlet case. To NLO, all splitting functions in ©-(0) have the perturbative 
expansion 



AP t n£,a s 



(u),(i) 



(0- 



(14) 



Just like their unpolarized counterparts, the polarized space-like and time-like splitting 
functions are equal in LO: 



Aif WO) (0=A^ MU; (0 



(S),(0), 



1.1 



1.1 



(15) 



Eqs. ([1~5|) are manifestations of the so-called Gribov-Lipatov relation (GLR) [18| which 
connects space-like and time-like structure functions within their respective physical re- 
gions (£ < 1) and is known to be broken beyond LO in the unpolarized case [15|, [II 
Recalling that for x < 1 @, 21 

1 + x 2 



where 



AP^\x) 
AP£>.<°>(*) 
AP^\x) 
AP^M ix) 



Cf — g, Ca — 3, Tf 



1 — x 

27} [2x - 1] 

C F [2 - x] , 
2C A r ' 



-2x + l 



x 



T R n f = -n f , fa = yC A - -T f , 



(16) 



(17) 



it becomes obvious that the space-like and time-like LO quantities are also directly related 
by analytic continuation through x — 1: 

1 



AP 



cn,(o), 



AP™\z) = -zAPg>>®£) , AP^(z) 



Cf_ 
2T f 



A p(S),(0)/iN 



AP^(z) 



22/ 



ApW.W(i), AP™°)(z) = -,Ap(f'(°)(i) 



'1? 



where z < 1. For future convenience we have explicitly written out the singlet as well 
as the non-singlet sector even though all LO quark-to-quark splitting functions coincide, 
AP^i (0) = AP g ( ^ {0) = APgpM. Eqs. (0) represent the analytic continuation or Drell- 
Levy-Yan relation (ACR) [[HJ to LO which we cast into the generic form 



AP<P 0) (z) = zAC 



ap}; 



(5),(0), 



X 



1, 

z' 



(19) 



where the operation AC analytically continues any function to a:— >l/z>l and correctly 
adjusts the color factor and the sign depending on the splitting function under considera- 
tion, cf. Eqs. fllSP . The LO relations ( fLg|) are based on symmetries of tree diagrams under 
crossing, and one therefore has to expect that they are in general no longer valid when 
going to NLO, depending on the regularization and the factorization/renormalization 
prescriptions used in the NLO calculation. This is exactly what happens in dimensional 
regularization in the MS scheme as was shown in |15|] for the unpolarized non-singlet case. 
Fortunately, as was also demonstrated in flI5| , the breaking of the ACR arising beyond 
LO is essentially due to kinematics and can therefore be rather straightforwardly detected 



within the method used in [15, 16, 22] to calculate splitting functions. We will now first 



collect the findings of [15[] concerning the connection between the space-like and time-like 
flavor non-singlet configurations in the unpolarized case and analyze in detail their ex- 
tension to the singlet sector made in |I6| . Afterwards we will apply the results to the 
polarized case. 
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3 NLO results for the unpolarized case 



Eqs. (f|)-([n]) above have been written down for the polarized case, but they all apply 
equally well to the unpolarized one when all A's are removed and, obviously, the unpo- 
larized LO splitting functions Pjp'^ as calculated in 



Cj 



1 + x 2 



1 — X 



Pg>>V>(x) 



27) 
C F 

2C A 



x 2 + {l 
1 + (1- 



x 



X) 



X 



X 



X 



+ x(l — x) 



(20) 



1 — X X 

(for x < 1), are used in Eq. (|T6|). Furthermore, in the unpolarized case with pure photon 
exchange one needs to introduce two independent structure functions J 7 ^, (see, 



e.g., p3[ ^4], ^5fl ) with short-distance cross sections C\ x , C, 



(U) MU) 



i,2 



q, g), respectively. 



In fll5| ) |16| , |22]| the unpolarized NLO evolution kernels for the space-like situation were 
calculated using a method that is as close as possible to parton model intuition since 
it is based explicitly on the factorization properties of mass singularities in the light-like 
axial gauge. The general strategy here consists of a rearrangement of the perturbative 
expansion which makes explicit the factorization into a part which does not contain any 
mass singularity and another one which contains all (and only) mass singularities. More 

explicitly, Mj^ (j = q, g, k = 1, 2), the contribution of virtual (space-like) photon-quark 

( s) 

or photon-gluon scattering to the structure functions T\ on parton-level, is expanded 
into two-particle irreducible (2PI) kernels. In the light-cone gauge these 2PI kernels have 
been proven to be finite as long as the external legs are kept unintegrated, such that 
all collinear singularities originate from the integrations over the momenta flowing in the 
lines connecting the various kernels [2(|. This allows for projecting out these singularities 
|T5| |, and Mj t k can thus be written in the factorized form 



M. 



(S) 



ij 5 



(21) 



1=1,9 



where the C\ S ^ are finite (and obviously depend on the hard process considered), whereas 



the r£p = rlp(x, a s , 1/e) contain just the mass singularities (which appear as poles in 



e when using dimensional regularization, d — 4 — 2e) and are process- independent . The 
are to be convoluted with bare ('unrenormalized') parton densities which must cancel 



their poles. As was shown in [15|, the MS scheme Altarelli-Parisi [£j kernels one is looking 



for, appear order by order as the residues of the 1/e poles in V^p , 



rjf(x, a s , -) = 5(1 - x)5ij - - 



<M p (5),(o) 1 '■" x 2 
27ry ij 



C (22) 



The NLO contribution to the hard short- distance cross sections in the MS scheme is 
obtained by calculating the full ('bare') subprocess cross sections C^P (x, -) (i = q,g, 
k = 1,2) and subtracting off the poles: 



c«r (*> = ^ (». i> + h^) " ^H. swo) m . (23) 

where /i is the arbitrary mass scale to be introduced in dimensional regularization. 



In the time-like region one can repeat the above procedure and introduce analogous 
quantities T^\z,a s , 1/e) that contain all final-state mass singularities arising in a frag- 
mentation process. It turns out [115] that the task of establishing the connection between 



(S) (T) 

T\j and via analytic continuation can be reduced to understanding the differences be- 
tween the 2PI kernels in the space-like and time-like situations. These essentially amount 
15fl to relative extra phase space factors of (k 2 • n/k\ • n)~ 2e in the time-like case, where 



k\ and k^ are the momenta of the particles entering or leaving a 2PI kernel, respectively. 
Here n is the vector specifying the light-cone gauge and the longitudinal direction, i.e., 
(k2-n/ki-n)~ 2e = (~ 2e with ( to be interpreted as the fraction of the momentum k\ trans- 
ferred to the particle with &2- In the unpolarized case a further difference arises from the 
spin-average factor for initial-state gluons which is id — 2) _1 = 1/2(1 — e) in d dimensions. 
As apparent from (0), the off-diagonal splitting functions interchange their roles during 
the transition from the space-like to the time-like situation. In particular, the space-like 
Pjp which includes the spin-averaging factor (d — 2) _1 gives rise to the time-like Pjjjp 
which should just have the spin- average 1/2, and vice versa for P[;p, Pg). These effects 
have to be taken into account along with those coming from the (k 2 • n/k\ ■ n)~ 2e terms 
mentioned above. Consequently, all this gives on aggregate for z < 1 

T (T) (7 rv -\- _ 7 l-2^r (S) (- rv -\ 
- L go,±^) u: S) ) * L qq,±\ ) u Si ) i 

e z e 



c V-*(i-e)rg>(i, 

2 e 



2T, 



2T f z 



l-2e 



C F 1-e 



rS } (-,« s ,-), rg)(z,o.,i) 



y l-2e F (5)/^ 



1 ^ 

z e 



(24) 



We also include now the corresponding relations for the hard subprocess cross sections 



Cii before subtraction of their pole terms (see Eq. (|23|)): 



_l-2eA(S),(l)/l 1 



c 



9,fc I 



1 1 



2T/ (1 e)C ^ fc 



(25) 



It becomes obvious that higher pole terms in the expression for rh- in ( p2j ) will generate 



additional contributions to the single pole of when they are combined with the factors 



-2c 



or (1 — e) 



±i 



in 



-2c 



\nz + 0{l) 



(26) 



In the same way the pole terms in c\ S ^ ,<yl \x, -) will give rise to extra finite contributions 

V(T),(1) 



to the QV that remain after the pole is subtracted. Following |L5| we separate all such 
ACR-violating contributions by writing 



T ( if\z,a s , -) 



zAC 
zAC 



-I h 

z e 



a 



i,k ' 



X 



1 1 

z e 



+ T e ij (z,a s ,^) , 



(27) 



where, as before, k = 1,2, i,j = q,g (or Hj = qq, ±' for the non-singlet case). We have 
extended the notation ,AC[...] for the analytic continuation (see Eq. (|19|)) to the short- 
distance cross sections, its action here being obvious from Eq. (^5j). One can now go 



through the NLO calculation [^] of the rj^'^ graph by graph to pick up the 1/e 2 pole 
terms and thus to extract the contributions to the T\Az, a s ,l/e) that break the ACR. For 
this purpose we present in Fig. 1 the basic topologies for all NLO diagrams involved here, 
where the notation is as introduced in |L5], For topologies (cd),(e),(fg) the higher 

pole terms are necessarily proportional to the pole terms in the renormalization constants 
as listed in |~5], EIJ. The corrections to the ACR coming from these graphs are therefore 



3 We note that the remaining topologies ((b), (jk)) introduced in |15|, £2| do not possess higher pole 
terms and thus do not contribute to the Tf-. 
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(h)-(l) (e) (cd) (fg) 



Figure 1: Basic topologies of the diagrams that contribute to rf,- (as defined in Eq. $2l\)) 
in NLO. 



easily detected upon insertion into Eq. (|23|). The higher pole terms of the genuine ladder 
graph (h), which is quadratic in the LO 2PI kernels, are given by the convolution of the 
LO splitting function for the upper rung with that corresponding to the lower one, and 
again the result is straightforwardly obtained. The only subtlety arises for the subtraction 
graphs (i) whose contributions to T^'^ are proportional to [|13[ 

graph(i) ~ i ((1 - zr4 S K) ® n ( f* (0) > (28) 

where -P^ 4I2I denotes the (d = 4 — 2e)-dimensional LO splitting function standing for the 
upper part of the diagram, the factor (1 — z)~ e arising from phase space. The lower part 
of graph (i) is represented by P^'^ which is the usual /bwr-dimensional LO splitting 
function. Application of the rules (E^) to topology (i) is then to be understood as to 
include the kinematical z~ 2e corrections only in the kernel representing the upper part of 



the diagram ||15|| , which gives 

'1 



graph (i) ~ zAC 



-I-21bz4 S, ' (0) )®Pt? M 



(29) 



as the contributions to the Y\y If required, the spin- averaging factors (1 — e) ±1 also have 
to be taken into account. Contrary to all other topologies, adjusting the spin-averaging 
in graphs (i) generates corrections also to the ACR for the diagonal NLO singlet splitting 
functions Pqqps an d Pffl'^. As an example, Fig. 2 shows the graph of topology (i) for 
the case of the 'pure singlet' function P^pp for both the space- like and the time-like 
situations. According to Eq. (p8"D , the contribution of this graph in the space-like case 
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pinm p(T),(o) 



Figure 2: Graphs of topology (i) contributing to P^ ps m the space-like (left) and time- 
like (right) situations. 

involves the convolution of the <i-dimensional splitting function Pifl'i ^ - which contains 
the spin-averaging factor (1 — e) _1 - with the /bt/r-dimensional function Pfjp'(°\ In the 

(T) (0) 

time-like case this obviously turns into the convolution of the rf-dimensional -P 9<?j 4_2 e with 
the /bur-dimensional P^^\ and no factor (1 — e) _1 is involved. The reversed situation 
appears in the CpTf part of P^'i 1 ). It is straightforward to account for these effects. 

We have now collected all ingredients for calculating the corrections to the ACR arising 
from the z~ 2e terms (and, if applicable, the spin-averaging factors) in Eqs. (|24]), i.e., for 
the functions rf- and Cf^ in (|27|) . Following Eqs. (|22|) , (|23|) we keep just the residues of 
the 1/e poles in rf . and subtract the pole terms from the 'bare' subprocess cross sections. 
We then rewrite Eq. (|2"T| ) as 



P$ m (z) 



C 



i.k I 



zAC 



zAC 



z 



a 



.k 1 



X 



1. 

Z' 



(30) 
(31) 



where, using the unpolarized LO splitting functions of Eq. (f20"|), we have 



P^2(z) = {3 P^(z)hiz , 



r qq,PS\ Z ) 



CfTj- 



9z 



[1 - z)(38 + A7z + 38z 2 ) + 3(1 + z)(4 + llz + 4z 2 ) lnz 



2{C F -C A )P^\z) 



+ 
+ 



1z 

CaCf 
9z 



45i(z) + 2 ln(l -z)\nz + ln(l - z) 
- 4z(3 - z) In z - 20 + 28z - 5z 2 

6 (4z 3 + 15z 2 + 18z + 19) lnz + (1 - z)(235 + 55z + 64z 2 
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-T f pW-V>{z){l-2\nz) 



+ 2{C F -C A )P& 



+ C F Tf 
2C A T f 



+ 



9z 



ASx(z) - 2 ln(l -z)\nz + ln(l - z 
-4(2 - z) (1 - 2z) In z - 5 + 28z - 20z 2 
6(-4 - z - 46z 2 + 9z 3 ) In z - 64 - 24z - 114z 2 + 169z 3 



99 



2\n zP^'W 



VI 



C;f(z) = (21nz + l)p(f(°)(z). 



(32) 



For the off-diagonal Ph in (R2I) we have introduced the function JT6 



with the Dilogarithm function Li 2 . 



i-» ln(l-y) 

= -Li 2 (l - -z) 



(33) 



The final step is to determine the analytic continuation of the space-like NLO splitting 



functions P^p'^ as published in [[15], 0, |22| and of the short- distance cross sections C, 



(see, e.g., p5| ) by using the operation ^4C[...] defined above. This is a straightforward 



task apart from two subtleties |T3]. Firstly, one has to recall that - as a result of the 

finiteness of the 2PI kernels in the light-cone gauge - the 1/e poles in the expression for 

(s) 

Tj { originate from the final integration over the momentum k of the off-shell particle 
emerging from the uppermost kernel (see Fig. 1 for notation). This momentum obviously 
satisfies k 2 < in the space-like case, but becomes time-like (k 2 > 0) when dealing with 

(T) 

the T\j . In other words, the transition from the space-like to the time-like situation does 
not only involve analytically continuing to x > 1 but also to k 2 > 0, crossing the threshold 
at k 2 = 0. Since the 1/e poles arise from terms like (— k 2 )^ 1 ^ 6 (for the space-like case) 
which are integrated over — k 2 down to —k 2 = 0, factors of (— l)~ e will appear in the 
virtual integrals when going to k 2 > 0. Using Re[(— l)~ e ] ~ 1 — e 2 7i 2 /2 one realizes that 
this will result in extra ^-contributions when multiplied by double pole terms present 
in intermediate stages of the calculations [|15||. The latter arise from the emission of 



soft gluons in the qqg vertex. A similar feature is present for the qq'j vertex p3] and 



thus appears in the quark short- distance cross section Cq k . Here it is the transition 



13 



q 2 < —>■ q 2 > that is responsible for the effect [0| which only affects the endpoint 
contributions at z = 1 to be discussed below. 

The other subtlety concerns the analytic continuation to x > 1 of terms ~ ln*(l — x) 
(i = 1,2) appearing in the Pjp (x) and in the C\ S k . To find the correct answer 
for this, one has to go through the relevant real (phase-space), virtual, and convolution 
integrals in the limit x > 1. It turns out (see also |15[) that all integrals can be smoothly 
continued through x = 1 via ln l (l — x) — > ln*(x — 1), the only relevant exception being one 
of the scalar three-point functions with a light-cone gauge propagator which, for x < 1, 



was given in Eq. (A. 15) of ref . [ 22| . This particular three-point function only contributes 



to the NLO splitting functions, but not to the hard subprocess cross sections C^^ . 
Upon recalculation of the function for x > 1 one finds that the correct continuation yields 
ln 2 (l — x) — > ln 2 (x — 1) + tc 2 in this case. When combining this result with that for the 
crossing of the k 2 = threshold discussed above, one arrives at the interesting finding 
that the correct analytic continuation to x > 1 of all terms ~ ln*(l — x) {i — 1, 2) in the 
NLO splitting functions PjP' ( " 1 \x) is effectively obtained by simply substituting 

In (1 - x) — ► In (x - 1) , (34) 
In 2 (1 - x) — ► In 2 (x - 1) - 7T 2 . (35) 

For the non-singlet case, in which no ln 2 (l — x) terms appear in the space-like NLO 



splitting function, this result is in agreement with the conclusion drawn in ref. |15| that 



the extra 7r 2 terms stemming from the threshold at k 2 = and from the three-point 
function cancel each other. In case of the NLO short-distance cross sections C^^ there 
are again only single powers of ln(l — x), and Eq. ( p4[ ) provides their correct analytic 
continuation. 

Combining everything, we arrive at the final result for z < 1 for the NLO (MS scheme) 

(T) (1) 

non-singlet and singlet time-like splitting functions (z) and the NLO time-like 

short- distance cross sections Cf^ (z). The result for the NLO splitting functions is in 
complete agreement with that of [|1^, |I6| , apart from a known misprint^ in [I6|. The 



4 Thc term (10 - l&c - 16x 2 /3)lna; in the C F T f part of P^ )A1) {x) in ^ must correctly read 
(-10- 18a; - 16a; 2 /3) In ir . 



14 



endpoint contributions, i.e., the terms ~ 5(1 — z), to the diagonal splitting functions can 
be obtained from the fermion number and energy-momentum-conservation conditions 
|15| , |l~6f and are exactly the same as in the space-like situation. In case of the NLO time- 
like quark short-distance cross section the endpoint contributions differ from those in the 
space-like situation by Cf7c 2 5(1 — z) which is just the effect of the above mentioned n 2 - 



correction when crossing the threshold at Q = [15]. Taking this into account, the results 



for the agree with those in, e.g., p5| . Since all the unpolarized expressions 



have appeared in the literature we do not repeat them here. We only note that, in 
contrast to the leading order (cf. Eq. (|TJ)), the NLO differences P;f ),(1) (0 - ^f ),(1) (0 

(T) (1) (S) (1) 

and C q J (£) — C q J (£) are non-zero (note that here £ < 1 in both the space-like and 
the time-like functions), i.e., in addition to the ACR the GLR is also broken beyond LO 
|IB|, as we anticipated in the introduction. 



4 The breaking of the ACR revisited 

Before addressing the polarized case which we are mainly interested in, let us return for a 
moment to Eqs. (|30|)-([32l). The rather simple structure of the P^ 1 \z) and its transparent 
origin suggest that there could be a more straightforward way of linking the time-like 
and the analytically continued space-like NLO splitting functions, than going through 
Fig. 1 graph by graph and picking up the higher pole terms. The starting point for such 
considerations is to notice that Eq. (^) (when adapted to the unpolarized case, i.e., with 
the A's omitted and the P,-f as given in fl20|)) only states that the /bur-dimensional 
LO splitting functions obey the ACR. The rule must break down for the (d = 4 — 2e)- 
dimensional counterparts, Pj^-^j °f ^ ne P^ as an immediate consequence of Eq. 
We write down a LO analogue of Eq. (RO) in d dimensions, 



P$%(z)=zAC 



r ji,i-2e\ Jj — , 

' Z 



6,(0), 

P e q f{z) = (~2\n z)eP^(z) 



+ P l f\z) (36) 



(for z < 1), the main difference being that the LO P^ (z) are only O(e) and not 0(1) 



5 There is a typographical error in the first equation of appendix II in |25[ : the prefactor of the In a; 
term should correctly read 3(l + a; 2 )/(l — x). 
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-2lnz + l)eP^ 



;-21n z)ePg>>®(z) 



(37) 



As already seen from the example of Eq. (|28|) , the pieces ~ e in the ci-dimensional LO 
splitting functions result in finite contributions in the calculation of the NLO splitting 
functions. One therefore suspects that the breakdown of the ACR beyond leading order 
in the MS scheme, as expressed by Eqs. (j30|) -(32), is entirely driven by the corresponding 
breaking in the part ~ e of the <i-dimensional LO splitting functions, cf. Eqs. (|36|) 

(T) (1) (T) (1) 

If this is indeed the case, then the functions P>, , C ik and 



c 



(T),(l) 



= zAC 
= zAC 



P^ m (x = -) 



z 

h 

z' 



(3* 



respectively, should be simply related by a factorization scheme transformation^. The 
general form of such a transformation can be determined from the condition that is must 
leave any physical quantity such as, e.g., or invariant, and reads 



P (n(i) 

p(T),(l) 



p 



^0 IT) 

2 qq ' 

p(T),(l) _ A>£(T) , 

2 



£(T) p(T),(0) 



a 



,k 



c 



(T),(l) 



ST) 
^iq 



(39) 



where the subscript '®' denotes convolution when evaluating the commutator. Again, 

p(T)M and p(T),(l) 

are the (unpolarized) LO and NLO evolution matrices, respectively 
(cf. Eq. (|T2D ) , and z^p and the 2x2 matrix generate the transformation. In analogy 
with Eq. (0) we setQ 



£(T)_ 



^qq 

2n f qg 



2 MV 

AT) 



(40) 



z 



yy 



6 For the non-singlet case this possibility was already hinted at in |l5| . 

7 For our purposes, we do not need to distinguish between the non-singlet and the qq entry in the 



singlet matrix Z^ T ) even though these could be chosen differently in principle. 
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According to Eq. ( p7|) one now expects that the choice 



(2lnz + a ij )P i 



(S).(0). 



'j 



(41) 



with the logarithms being of kinematical origin and the resulting from the adjustment 
of the spin-averaging factors, 



>-qq 



a 99 ~ ' 



1 



-1 



(42) 



transforms all time-like NLO (MS) quantities to a scheme in which they satisfy the ACR, 
i.e., in which 



P, 



r (T),(l) 



P, 



a 



i,k 



zAC 
zAC 



C 



(S),(l) 
i,k 



X 



1. 



(43) 



This indeed turns out to be the case as one finds upon insertion of the z\a in (|41|) into 



Eq. (|3^) and comparison with (0). We emphasize that the space-like NLO quantities on 
the right-hand-sides of Eq. ( f43l) have not been transformed and are still in the MS scheme. 
Eq. fl4*3| ) therefore links quantities referring to different factorization schemes. This is per- 
fectly legitimate since one is free to choose the factorization schemes independently for the 
space-like and time-like casesQ. On the other hand, it does not really appear sensible from 
a physical point of view, and it actually turns out J^J that the transformed time-like 
NLO splitting functions of fl43| ) do no longer obey the energy-momentum-conservation 
condition. Anyway the above scheme transformation is not meant to be used in any 
practical calculation, it just serves to identify the breakdown of the ACR beyond LO as 
a mere matter of convention and provides a very transparent and remarkably simple way 
of obtaining the correct MS time-like splitting functions from the analytically continued 
space-like ones. We note that in |H| the unpolarized NLO time-like splitting functions 
and short- distance cross sections were calculated using the cut vertex method. In this 
formalism the validity of the ACR occurs quite naturally if certain renormalization condi- 
tions are imposed [pq] , and the results of [19j therefore correspond to the P^ , C ik 
in fl43|) rather than to the MS scheme results. 

8 For instance, one can choose to factorize initial- and final-state collinear singularities differently in 
any higher order calculation. 
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5 NLO results for the polarized case 



The extension of our results in sections 3 and 4 to the spin-dependent case is rather 
straightforward now. We first write down Eqs. (|30|) , (|3T|) for the polarized case, 



AC. 



cn.(i) 



zAC 
zAC 



AP^\x 



AC 



(S).(l), 



X 



+ AP^\z) 



+ AC 



(44) 
(45) 



where z < 1 and where we have recalled from section 2 that contrary to the unpolarized 
case there is only one longitudinally polarized structure function for pure photon exchange, 
G\ ■ For the space-like situation, NLO (MS) results for the short- distance cross sections 



AC} 



, ( (i.e., the coefficient functions for g\ ) have first been published quite some time 
ago |29|, [3i~l , whereas the corresponding MS splitting functions have been calculated 
only fairly recently via the OPE |L2] (where they appear as the anomalous dimensions) 
141, where the method of 



and in 



was used. To be more precise, use of 



dimensional regularization in such NLO calculations for the polarized case implies to 
choose a prescription for dealing with the Dirac matrix 75 and the Levi-Civita tensor 
e^ vpcr in d 7^ 4 dimensions, which poses a non-trivial problem. In |12|] the 'reading point' 
prescription of |£J with a fully anticommuting 75 was chosen, whereas JT4J adopted the 



original definition for 75 of |33| (HVBM scheme) which is widely considered to be the most 
consistent method. A crucial feature in both |12| and [13|, [L4| was that the genuine ('naive') 
MS scheme result for the non-sine let NLO splitting function AP^ (1 \x) (cf. Eq. (0)) 



possessed the disagreeable property of having a non-zero first moment (x-integral) , in 
obvious conflict with the conservation of the non-singlet axial current f34] , |35| which 

(S) 

demands that the first moment of the non-singlet quark combination Aq\_ be independent 
of Q 2 . This effect was clearly due to the 75 prescriptions chosen and could be removed by 
a finite renormalization in |12j or, equivalently, by a factorization scheme transformation 
in [13, 14] generated by the difference of the d-dimensional LO quark-to-quark splitting 
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functions for the unpolarized and polarized (HVBM scheme) casesfl, 



Uir qq,A-2e y 



X)-P t 



<?<?,4-2e 



X) 



AC F e(l - X) . 



(46) 



It turned out that both |12[ and |U| [L4[] then arrived at the same final result for the space- 
like polarized NLO splitting functions and also for the coefficient functions ACf^'^ 1 ' and 
ACf for which the previous results of [E5|, [JU] and |JT], respectively, were confirmed. 



In a strictly technical sense of the word, the results of |12], [L3[ [L4[] thus do not correspond 
to the MS scheme. On the other hand, the '75-effect' described above has been known 

39| and is purely artificial 



35 



to occur in the HVBM scheme for some time |37|, [ 
in the sense that it is related to helicity non-conservation at the quark-gluon tree-level 
vertex in d 7^ 4 dimensions as expressed by the non- vanishing of the rhs of Eq. ([46]). 
Since physical requirements like the conservation of the non-singlet axial current serve to 

are 



remove the effect in a straightforward and obvious way, the final results of [|12|, [13], |I4 
nevertheless usually regarded as the 'real' conventional MS scheme results. 



The reason for going into this discussion is the following. If we use the final results of 
fl2| , [13| , |l~4"|| for the APjp'^ and the AC^ to obtain their time-like counterparts via 



Eqs. (|44]) , (|45|) , the factorization scheme transformation generated by ( f46|) and performed 
in the space-like situation will obviously also affect the time-like result. On the other hand, 
in the case of the time-like NLO quantities, there appears to be no obvious requirement 

1 (T) (1) 

that enforces a certain value for, say, the first moment J APg q + (z)dz. Thus in principle 
one would be allowed equally well to use the space-like 'naive' MS scheme results in 
(fB]),(fI^), i.e., those that possess the wrong (non- vanishing) value for the integral of 
APlf ]^\x). However, taking into consideration the origin of the above '75-effect' as a 



pure artefact of the dimensional calculation in a certain 75 prescription, we decide against 



this latter option and will use the final results of [12, O, 14|, i.e., the 'real' MS scheme 



results for the space-like case in what follows. This choice of factorization scheme appears 
most sensible since it actually turns out that Eq. (K6|) remains completely unchanged 



9 As was also shown in [|l3| [l4| p5| , the scheme transformation corresponding to ( puf ) is needed at the 
same time to bring the first moment of the quark non-singlet coefficient function into agreement with the 
value given by the Bj0rken sum-rule pa. 
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r qq,i-2e\ z ) 



when transformed to the time-like situation, 

4CVe(l - z) , (47) 

implying that the unphysical helicity non-conservation at the quark-gluon vertex in d 7^ 4 
dimensions also takes place in the time-like case if one uses the HVBM prescription for 75 



TTH . Our choice obviously has implications for factorizing collinear singularities in NLO 



calculations of other cross sections with polarized final state particles: The 'real' MS 
scheme factorization counterterm for all collinear poles coming from polarized (time-like) 
quark-to-quark transitions should be taken as (llj 

where Mp is the factorization scale and Acrf^ 2(E is some appropriate polarized Born-level 
cross section in d dimensions. 



After these precautions we can turn to the analytic continuation to x > 1 of the 
space-like NLO quantities which is required for Eqs. ([44]), ([45]) and works in exactly the 
same way as for the unpolarized case studied in section 3. The other ingredients to 



Eqs. (^4|) ,d45|), AP^ (z) and AC* (z), are also straightforwardly calculated following 
the lines of section 3. The situation is facilitated by the fact that in the polarized case 
there are obviously no complications due to the gluon spin-averaging factors. Thus we 
only have to keep track of the z~ 2e terms. We find: 



>«,(!) 



AP^\z) 
ACf\z) 
AC^(z) 



,6,(1), 



(3 AP^°\z) hxz = P e qq ^ 
-12C F T f [(1 + z)lnz + 2(l - z)} , 
4(C F - C A )ApW><®(z) [-2S 1 (z) + ln(l - z) hxz] 



+ C 



(4-5z) In z -4(1 



8C A Cf 



[1 + z) \nz + 2(l-z) 



T f ApVp>W(z) \nz + 4(CV - C A )AP^°\z) [2^ (z) - ln(l - z) \nz 



+ 2C F T 



(5 - Az) \nz + 4(1 - z) 



+ \c A T, 



lOz-23) In z- 24(1 -z) 



AP; q %(z)+(3 AP^°\z)lnz 



2\nzAp(p' {0 \z) 
2\nzAP^°\z) 



(4£ 
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As for the unpolarized case in section 4 it turns out that all these terms can also be 
fully accounted for by a factorization scheme transformation, i.e., the MS scheme 

(T) (I) 

and AQ and the corresponding analytically continued space-like functions, 



A p(T),(i) 
AC- 



= zAC 
= zAC 



ap; 



(S),(l) 



AC 



(S),(l) 



1* = -; 

(» = i; 



(49) 



respectively, are related via Eqs. ([39|),([40|) (with A's everywhere in (|39|),([40|)) if one chooses 



21nzAP 



(S),(0), 



'j 



(50) 



Now we finally insert everything into Eqs. (|44|) , (|45|) and arrive at the final results for 
the time-like NLO quantities which, in case of the splitting functions, are conveniently ex- 
pressed as differences with respect to the space-like situation, at the same time indicating 
the breakdown of the GLR: 



AP 



(T),(l) 



APf ^(z) + A t ,(z) 



(51) 



where the AP- 



Aqq,Ps(z) 



A 99 (^) 



z) are found in [ I2"| , [13], |I4] and 
In z 



z 



l + 3z 2 ) lnz + z 2 + 4z + 1 



4 (l + z 2 ) ln(l - z) 
4 C F T f (In z - 3) [(1 + z) In z + 2(1 - z)} , 

(2 - z) f-245i(z) - 4tt 2 + 6 ln 2 (l - z) + 12 ln(l - z) In z 



(52) 
(53) 



~Cf 
3 F 



-31n 2 z + 15 In z 



+ (15z - 6) ln(l - 2) - 33z + 54 



+ - C A C F [6(2 - z) [l2Si(z) + 2vr 2 - 3 ln 2 (l - z) 
9 l 

+ (6 - 39z) ln(l - 2) - 18(^ + 4) In 2 z + 72(3z - 1) In z - 71z + 4 

+ ^C F Tf [3(2- z)ln(l- z) + z + 4] , 
9 
g 

= - Tf [—3(22; — l)(ln(l — z) + lnz) - 4z - 1] 
9 

2 r r 

+ -C F T f (2z-l) 24Si(z) +4tt 2 - 61n 2 (l - z) - 3 In 2 2 
3 

+ (6z - 15) ln(l - z) + 30 In z - 78z + 57] 



(54) 



+ -C A T f [6(2z-l) 



-12£i(z) -2vr 2 + 61n(l-z) lnz 
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+3 ln 2 (l - z)] - (Qz - 39) ln(l - z) + 36(1 + z) In 2 z 
-3(26^ + 11) Inz + 284z - 217] , 

= AC F T f [6(1- z) + 6\nz + (1 + z)\n 2 z 

8 2^2-3^+2 
3 ' 1-z 



+ ^1 



In z 



3 "l-z 
-26z 2 + 63z - 26 



For the short-distance cross sections we obtain 

'ln(l -z)\ 



Yl(2z 2 - 3z + 2) ln(l - z) + 6(3z - 4) \nz 



ACjr^\z) = C F 
+ 



1+V 



1- z 



+ 2- Inz — r 



l-z 



2(1-*)- 



(2-z) ln(z 2 (l - z)) -4 + 3^ 



in agreement with the results of [|TTJ for the corresponding choices AfF 1 



(55) 



(56) 



(57) 
(58) 

-4(1 - z) 



and Af®(z) = in Eq. (14) of that paperQ. The '+ '-prescription in Q57| ) is defined as 
usual via 

1 dzf(z) (g(z)) + = f dz (f(z) - /(l)) g{z) . (59) 
o Jo 



For numerical evaluations of our results it is convenient to have the Mellin-moments 
of the expressions above which are defined by 



f[n) 







f(z)dz 



(60) 



and are presented in the appendix. Fig. 3 provides a comparison of our results for the 
NLO (MS scheme) time-like polarized and unpolarized splitting functions in Mellin-n 
space. It is interesting to observe that all LO and NLO time-like splitting functions obey 



AP™*>[n] 



P-j [n] (k = 0, 1; i, j = q, g) as n — > oo, i.e., as z — > 1, except for 
^■Pgq [ n ]- A similar observation was made for the space-like quantities where, again, 
only APjf - (1) [n] does not fulfil AP l 



f -«[n] - ' (fc) [n] asn 



oo g|. Finally, the 



10 Note that our definition for the gluonic short-distance cross section differs by a factor of 2 from the 
one used in 1 1 . 



22 



10 

-10 
-20 
-30 

20 
10 

-10 



- / 



-_ I 

; / 
- I 



n 1 1 — i — i — n~ 



t( ' (A)P^' (1) [n] 
1 1 — I I I I I I 



(A)Pj W1) [n] 




_i i i i i i i 



10 



n 




10 



n 



Figure 3: Comparison of the spin-dependent NLO (MS) time-like singlet splitting func- 



(T) (1) 

tions AP/j [n] as functions of Mellin-n according to Eqs. (A.1)-(A.6) with the corre- 
sponding unpolarized ones, as taken from [p^j, for / = 3 flavors. 



values for the first (n — 1) moments of the polarized NLO singlet quantities are given by 



AP^'W [1] = 3C F T f , AP^)' W [1] = -3C A C F + ir 
AP$P>V[1] = 2T f p , AP(J)-W[1] = -\ C \ - AC F T f - \c A T f - ^C A p , 
ACf =0 , ACP>W[1] = -^C F . (61) 



6 A supersymmetric property of the NLO time-like 
splitting functions 



In this section we finally very briefly address a relation that is conjectured for an M = 1 
supersymmetric Yang-Mills theory and connects all singlet splitting functions in a remark- 
ably simple way in the limit C F = Nc = 2Tf = N (cf. ||41|| ). In, e.g., the unpolarized case 
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it reads 



6[pW)M(t)] = P g f>«(£) + Pj^«(£) - P£°' W (0 - P% U H® = • (62) 

In LO (i = 0), the relation is satisfied by the both the unpolarized and the polarized 
splitting functions and, trivially, in both the space-like and the time-like situations. Be- 
yond LO, one can expect it to continue to hold only if the regularization method adopted 
respects supersymmetry. One therefore anticipates that the NLO (MS) space-like and 
time-like splitting functions of dimensional regularization will not satisfy fl62~D, in agree- 



ment with the findings in [|16| [T2| . However, one regularization method that is applicable 



to supersymmetry is dimensional reduction ||2], a variant of dimensional regularization. 
The scheme essentially consists of performing the Dirac-algebra in four dimensions and 
of continuing only momenta to d (d < 4) dimensions. In order to match the ultraviolet 
(UV) sectors of dimensional regularization and dimensional reduction, specific countert- 
erms need to be introduced [^, |44| in the latter which include a finite renormalization 
of the strong charge. Once this is done, all remaining differences between the results 
for a NLO quantity in dimensional regularization and in dimensional reduction can only 
be due to the effects of mass singularities. They are fully accounted for [03, 



by 



the differences between the <i-dimensional LO splitting functions (as to be obtained in 
dimensional regularization) and the /bttr-dimensional ones (corresponding to dimensional 
reduction) . In other words, the breakdown of the supersymmetric relation for dimensional 
regularization is entirely blamed on the breakdown of this relation in the ~ e-parts of the 
rf-dimensional LO splitting functions of dimensional regularization. This feature was ex- 
ploited in [L4[] to transform the space-like unpolarized and polarized NLO splitting 
functions of MS dimensional regularization to dimensional reduction via a simple factor- 
ization scheme transformation and to establish the validity of Eq. ( |B2] ) for the obtained 
quantities^]. We will now extend the considerations of f45[ , T%\ to the time- like situation. 
In the unpolarized case, one finds for the NLO MS splitting functions of dimensional 
regularization in the limit Cp = Nc 

S[pW).W( z )] =N 2 



2T f = NM 



2 13 5 1 

1 h -z - z 2 + (-1 + 2z + Az 2 ) \nz 5(1 

3z 6 3 v ; 2 V 



z 



(63) 



11 The supersymmetry relation for the space-like NLO kernels was proved prior to |4^, |Tj] in the OPE 
calculations of Jl6[ [l2| in which the transition from dimensional regularization to dimensional reduction 
occurs as a finite renormalization. 
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where we have included the endpoint contribution. The corresponding expression for the 
polarized NLO splitting functions can be obtained from Eqs. (pl|)-([56|) and the space-like 
results of §0 0, 



11 7 1 
---z+(l-z)\nz--S(l 



(64) 



What we need to do is to perform factorization scheme transformations fl39|) of these re- 
suits, with the (A) z\a to be determined from the parts ~ e of the time- like <i-dimensional 
LO splitting functions as obtained in dimensional regularization. For our purposes we 
only need to consider (|39|) in the combination appearing on the lhs of Eq. ( |62j) and in the 
super symmetric limit: 

5[( A)p(n(i)] _^ 5[(A)P^]-^5[(A)Z^] + ((A)Pg^ + (A)P^)®5[(A)Z^] , 

(65) 

where now (3q = 3N. The calculation of the parts ~ e in the d- dimensional LO time-like 
splitting functions yields^ 

5[Z {T \z)} = N (l-2z + 2z 2 -±5(l-z)j , (66) 

5[AZ {T \z)} = n(i-z-^8(1-z)^ . (67) 

Upon insertion of Eqs. (|63|) , (|66|) (or (|64]),(|67D) into ( |65D one finds that the resulting trans- 
formed 5[P( T )'W] (or 5[Ap( T )'^], respectively) vanishes identically. We thus have demon- 
strated the validity of the Af = 1 supersymmetric relation also for the time-like NLO 
polarized and unpolarized evolution kernels in dimensional reduction. Apart from being 
interesting of its own, this finding also provides evidence for the correctness of our results 
in sections 3-5. 



7 Summary 



We have presented a calculation of the unpolarized and polarized time-like NLO splitting 

functions, needed for the NLO Q 2 -evolution of (spin-dependent) fragmentation functions. 

12 As discussed in the previous section the difference ( |47| ) between the LO polarized and unpolarized 
time-like quark-to-quark splitting functions arising in the HVBM scheme is already accounted for in 
Eqs. ©-(H). To obtain the result in @ one therefore has to use AP^[%(z) = P^^z) (see also 

0)- 
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The starting point for our considerations were |15|, [H], [13], [14| in which the corre- 
sponding space-like quantities were calculated within a method based on the factorization 
properties of mass singularities in the light-cone gauge. As was shown in |15| for the 
non-singlet case one can then determine the time-like counterparts via analytic continu- 
ation to x > 1, which is also the way we have pursued. It turned out that beyond LO 
there are certain terms arising from phase space (and, for the unpolarized case, from the 
gluon spin- averaging in d ^ 4 dimensions) which prevent the analytic continuation rela- 
tion (ACR) of |17| between the space-like and time-like splitting functions from remaining 
intact. The same statement applies to the connection between the space-like and time- 
like short distance cross sections of electroproduction and e + e~ annihilation, respectively. 
Nevertheless, the corrections to the ACR are rather straightforwardly calculable within 



the method of [I5|. Even more, we were able to show that in both the unpolarized and the 
polarized cases one can transform the results to a factorization scheme, different from the 
MS scheme, in which the breakdown of the ACR does not occur. In the unpolarized case 
our final MS results confirm those of ]15], [16| obtained within the same method, whereas 
in the polarized case our results are entirely new. Finally we have shown that, when 
transformed to dimensional reduction, both our unpolarized and polarized results for the 
time-like NLO splitting functions satisfy a simple relation motivated from supersymmetry. 
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Appendix 

In this appendix we present the Mellin-n moments (as defined in Eq. (|60D ) of our NLO 
results for the polarized case in Eqs. The corresponding n moments for the 

unpolarized time-like splitting functions and short- distance cross sections can be found in 
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|27|. As in flSip we write 



Aif Ml) M 



Ay[n] 



(A.l) 



where the Ai^f [ n] can be found in a compact analytic form0 in ||0| and the A.y[ 
are given by: 

A OT [n] = Cr, 2- 



C£ 45! (n 
4(7^1) 



3n 2 + 3n + 2 
n(n + 1) 
n + 2)(3n + 1) 



. . 7T 2 2rz + 1 
-25 2 (n) + — + 

3 n^{n + l) 2 



A 9g[ n } 



n 3 (n + l) 3 
n + 2 / 



(A.2) 
(A.3) 



+ 



n(n + 1) \ 
7n 5 + 22n 4 + 7n 3 



5f(n) + 5 2 (n) + 



7T 



3n 3 + n 2 - 18n - 8 
n 2 (n + l) 2 



«Si(n) 



24n 2 



22n-4 



+ CaCf 



n 3 (n + l) 3 
" + 2 , „ 2 



n(n + 1) 



ww-3&( B )) + (n t.'y, i "T l) g,w 



3n(n + l) 1 



(67n 4 + 10ln 3 + 34n 2 + 144n + 72) (n + 2) 
9n 3 (n + l) 3 
n + 2 „ , % (n + 2)(5n + 2) 



9 



n(n + 1) 



«Si(ri 



° ^2 

9 ' 



+ 2C F T 



n 



n(n + 1) 
n — 1 



Si(n) - 



n(n + l) 2 
5n 3 - 3n 2 + 7n + 3 
n 2 (n + l) 2 



.f 



n(n + 1) 
7n 5 + 7n 4 - 5n 3 + 5n 2 + 4n - 2 



3« 2 -I- ^ 

Sl(n) + 3S 2 (n)) + , _ xo «Sx(n) 



n(n + l) 2 



n 3 (n + l) 3 



+ - g C A T f 
lln 3 



n — 1 



n(n + 1) 
12n 2 + 37n + 12 



(3S 2 (V) -35 2 (n) -7T 2 



5i(n) 



+ 



n 2 (n + l) 2 
67n 5 + 34n 4 + 32n 3 + 98n 2 + 249n + 72 



n 3 (n + l) 3 



n 3 + 3n 2 



1 + t °AT f 



F f n 3 (n+l) 3 ' 3 

lln 2 + lln + 24 



AS 1 (n) 



3n(n + 1) 



„ , , 7r 2 2n + 1 

-25 2 n +- + 4 

3 + 

2n + 1 



7T _ 



-2<S 2 (n) + — + 4 



n 2 (n + l) 2 



(A.4) 



(A.5) 



(A.6) 



13 Note that the results presented in [M need to be divided by —8 in order to bring them to our 
normalization for the NLO splitting functions. Also note that the definition of the i-components in the 
non-sing let sector (see Eqs. (§-©) occurs in a reversed notation in jiOj, i.e., as =F. 
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with 



n i 

3=1 J 



(A.7) 



The analytic continuation of the Sj., required for a numerical Mellin inversion, is well- 
known [[|7j]. For the moments of the polarized NLO time-like short-distance cross sections 
we find 



ACf)>«[n] 



Sl{n) + hS 2 {n) 



1 



2 n(n+l) 



1 2 9 



n rr 



(n + 1) 2 2(n+l) 
(A.8) 



ACf)>«[n] 



2+n . , 4 4 3 3 
+ — — + 



n(n + 1) 



n n + 1 (n + 1 N 



(A.9) 
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